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We study why gold forms planar and cage-like clusters while copper and silver do not. We 
use density functional theory and norm-conserving pseudo-potentials with and without a scalar 
relativistic component. For the exchange-correlation (xc) functional we use both the generalized 
gradient (GGA) and the local density (LDA) approximations. We find that planar Au n structures, 
with up to n = 11, have lower energy than the three-dimensional isomers only with scalar- relativistic 
pseudopotentials and the GGA. In all other calculations, with more than 6 or 7 noble metal atoms, 
we obtain three dimensional structures. However, as a general trend we find that planar structures 
are more favorable with GGA than with LDA. In the total energy balance, kinetic energy favors 
planar and cage structures, while xc-energy favors 3D structures. As a second step, we construct 
cluster structures having only surface atoms with Oh, T<j, and L symmetry. Then, assuming one 
valence electron per atom, we select those with 2(1 + l) 2 electrons (with I integer), which correspond 
to the filling of a spherical electronic shell formed by node-less one electron wave functions. Using 
scalar relativistic GGA molecular dynamics at T = 600K, we show that the cage-like structures of 
neutral AU32, Auso, and Aui62 are meta-stable. Finally, we calculate the static polarizability of the 
two lowest energy isomers of Au„ clusters as a means to discriminate isomers with planar (or cage- 
like) geometry from those with compact structures. We also fit our data to a semi-empirical relation 
for the size dependent polarizability which involves the effective valence and the kinetic energy 
components for the homogeneous and inhomogeneous electron density. Analyzing that fit, we find 
that the dipole polarizability of gold clusters with planar and cage-like structures corresponds to 
the linear response of 1.56 delocalized valence electrons, suggesting a strong screening of the valence 
interactions due to the d-electrons. 

PACS numbers: 36.40.Cg, 36.40.Qv 



I. INTRODUCTION 

Small clusters of metal atoms behave differently than 
the bulk matter, because each additional atom, or even 
each additional electron, can drastically change their 
electronic and geometrical properties 1 -. The noble metals 
clusters, with valence electron filling nd 10 (n+ l)s 1 , differ 
from the simple s-orbital alkali metals^, but also present 
striking differences among Cu, Ag, and Au 4 ^. Well es- 
tablished structural differences among Au^ [y = 0,±1) 
and Ag^ or Cu^ clusters are the following: i) Au^ clus- 
ters, specially the anions (v — —1), adopt planar struc- 
tures up to larger sizes than Ag and Cu clusters, as 
demonstrated by combined experimental and theoretical 
studiea^iL&S; ii) experimental photo-electron spectra for 
noble metal clusters with 55 atoms^S indicate that sil- 
ver and copper adopt some symmetry, preferably icosa- 
hedral, whereas the pattern for AU55 corresponds to an 
amorphous structure 1 ' 1 '; Hi) anionic and neutral AU20 
show a tetrahedral Td geometr y 4 ! 12 ; 13 , but Ag 2 o and 
CU20 have amorphous-like compact C s structures; iv) an 
icosahedral cage-like structure has been found to be very 
stable for Au32 1 ^^, but not for silver and coppeni 4 ^. 



Other meta-stable cage-like structures for gold clusters 
have been proposed recentljiiSiii. 

The differences between Au and other noble metal clus- 
ters are usually attributed to relativistic effects^, which 
stabilize the 6s orbital and destabilize the 5d one, fa- 
voring the hybridization of these orbitals. However, al- 
though Pt shows as strong relativistic effects as Au^, it 
has been shown that the competition between planar and 
3D structures of Pt clusters is not affected by relativistic 
effects^. Notice that Pt 7 is three-dimensional (3D)i2*2fi, 
like Ag7 or CU7, but Au 7 is planar. Notice also that the 
largest s-orbital contraction due to relativistic effects oc- 
curs in Au^i. Consideration of the spin-orbit coupling 
does not alter the relative stability of scalar relativis- 
tic structures of Au„ clusters with n < 20, but it in- 
creases the binding energy by about 0.08 eV/atom (1.85 
kcal/mol)22. 

Comparison of density functional theory (DFT) results 
for Aug and Aus at several levels of theory (that is, 
different exchange-correlation functionals), with results 
from quantum chemical calculations using second-order 
perturbation theory (MP2) or coupled cluster methods 
(CCSD(T)), indicate that DFT predicts planar struc- 
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tures, but MP2 and CCSD(T) predict the lowest energy 
Aug isomer to be non-planar by 26.6 kcal/mol and 1.5 
kcal/mol, respectively^. Another recent calculation 2 ^, 
using ab-initio correlated-level theory, predicts Aus to 
be planar. 

Concerning DFT calculations, we notice that the type 
of exchange-correlation (xc) functional has a decisive 
influence on the structural properties of gold clusters, 
but it is not so critical for silver and copper clusters 2 ^. 
Thus, first principles calculations by means of the Siesta 
codeaS with scalar relativistic pseudopotentials, found 
planar structures of neutral Au„ clusters with up to 6 
atoms using the local density approximation (LDA) for 
the xc-functionalS, and up to 10 atoms using the gen- 
eralized gradient approximation (GGA)^. Such a dis- 
imilar result was corroborated by a variety of different 
pseudopotential and all-electron scalar relativistic DFT 
calculations using LDA 28 i 29 and/or GGA 29 i 30 . 

In a recent paper, Gronbeck and Broqvist 31 compared 
the different contributions to the binding energy of sev- 
eral planar and 3D structures of Aus and Cus clusters, 
optimized within GGA and LDA. They found that planar 
Aus isomers have a significative smaller kinetic energy 
than 3D ones, which was attributed to d-electron der- 
ealization. A correlation between strong s-d hybridiza- 
tion and high stability of planar structures was found in 
ref. |a for noble metal heptamers, but does not appear 
to be a general tendency of small Au clusters. Instead, 
the preference of planar configurations for Aus isomers 
was attributed to a sizeable d-d overlap and to d-electron 
delocalization^i. 

Recently, the cr-aromaticity in saturated inorganic 
rings was examined 3 —. Evidence for d-orbital aromaticity 
in square planar noble metal clusters^ 3 - and in triangular 
gold rings 34 was presented also recently. The spherical ir- 
aromaticity of I/j symmetrical gold fullcrcncs fulfilling a 
generalized 2(7+ 1) 2 s-electron rule 3 ^, led Johansson 1 ^ to 
explain the extra stability of a cage- like AU32 cluster with 
Ih symmetry compared to other space-filling isomers. 

A rough estimation of the photoabsorption response 
of several isomers of AU32 and Au42^&, suggests that the 
cage-like structures could be clearly distinguished from 
space filling isomers in optical absorption experiments. 
However, it is difficult to separate in these spectra the 
spectral features due to symmetry from the features due 
to empty-cage effects. In this paper we will compare the 
calculated static dipole polarizability of planar and cage- 
like Au n clusters with those for 3D and compact isomers. 

In section [H] we outline the first principles method 
used in our calculations. Results are presented and dis- 
cussed in section IIIII In subsection IIII Al we compare 
GGA and LDA equilibrium structures of Au„ clusters 
with 6 < n < 9. We will focus on the relation of 2D or 
3D isomers with the derealization of d electrons, follow- 
ing the ideas of Gronbeck and Broqvist 31 . In subsection 
IIII Bl we contrast again LDA and GGA predictions for 
the stability of magic cage- like structures Auig, AU20, 
AU32, AU50, and Aui6o, compared to amorphous- like fill- 



ing space isomers. The stability of these cage-like struc- 
tures against molecular dynamics at constant tempera- 
ture (600 K), and against the loss or gain of one electron 
is also tested. In subsection IIII CI we investigate the use 
of the calculated static polarizability as a physical prop- 
erty sensible to the cluster structure. In section I1VI we 
will present our conclusions. 

II. COMPUTATIONAL PROCEDURE 

We use the first-principles code Siesta2& to solve fully 
self-consistently the standard Kohn-Sham equations^ of 
DFT within the GGA as parametrized by Perdew, Burke 
and Ernzerhof 38 , and within the LDA as parametrized by 
Perdew and Zunger— . For each xc-approximation we use 
a norm conserving scalar relativistic pseudopotential 4 ^ in 
its fully nonlocal formii, generated from the Au atomic 
valence configurations 5d w 6s 1 6p° , and core radii which 
we have tested and reported in previous work s 4 ' 42 . Flexi- 
ble linear combinations of numerical (pseudo) atomic or- 
bitals are used as the basis set, allowing for multiple-^ 
and polarization orbitals. In order to limit the range of 
the basis pseudoatomic orbitals (PAO), they are slightly 
excited by a common energy shift (0.01 eV in this work), 
and truncated at the resulting radial node^ 3 .. In the 
present calculations we used a double-^ 5p, 6s-basis, with 
maximum cutoff radius 7.62 Bohr. The basis functions 
and the electron density are projected onto a uniform real 
space grid in order to calculate the Hartree and exchange- 
correlation potentials and matrix elements. The grid fine- 
ness is controlled by the energy cutoff of the plane waves 
that can be represented in it without aliasing (120 Ry in 
this work). 

To obtain the equilibrium geometries, an uncon- 
strained conjugate-gradient structural relaxation using 
the DFT forces^ was performed for several initial clus- 
ter structures (typically more than ten), suggested by 
the several geometries for Au„, Au~, and Au+ isomers 
obtained previously^. 

The static dipole polarizability of a cluster can be ob- 
tained by using the standard numerical finite field per- 
turbation method, in which the field dependent energy 
is expanded with respect to an external uniform electric 
field F, 

E = E°- mFi - ^aijFiFj (1) 

where i,j are cartesian coordinates and the dipole mo- 
ment and the static dipole polarizability are obtained 
as energy derivatives, fa = — ^-\f=o, an d cty = 

— of of l F =Q' respectively. The external electric field val- 
ues used in our calculations were (in a.u.) |F| = 0.000, 
0.001, 0.006, 0.010, 0.014, and 0.018. The energies calcu- 
lated for these values were fitted to a polynomial expan- 
sion to obtain the first-order and second-order derivatives 
of energies with respect to the electric field strength. The 
mean polarizability is calculated as, a — Tr(ay)/3. 
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FIG. 1: (Color online) Equilibrium geometry of the lowest 
energy isomers of gold clusters having planar (Au n , upper 
row) or three dimensional (Au^, lower row) geometry as re- 
sulting from LDA and GGA scalar relativistic calculations. 
For GGA the ground state is the planar Au n isomer. For 
LDA the ground state is the 3D Aujj isomer, except for n — 
6, whose structure is also the planar Au6. 



III. 



RESULTS AND DISCUSSIONS 



A. Planarity and d-electron derealization in Au„ 

The onset of three-dimensional structures of neutral 
Au„ clusters was calculated at n — 11 within GGA—, and 
at n — 7 within LDA 27 . In this work we compare the pla- 
nar and 3D lowest energy isomers of Au„ (6 < n < 9) cal- 
culated within the Siesta code2&, and using LDASi and 
GGA 38 for xc-functionals (with the corresponding LDA 
and GGA scalar relativistic pseudo-potentials). Figure 
shows our results for the geometry of these isomers. The 
lower energy isomer is the planar one for GGA, and the 
3D one for LDA, except for n — 6. Table [I] gives the 
binding energy difference between 2D and 3D isomers, 
AE b = E b (Au n ) — E b (Au* l ), for the relativistic and non- 
relativistic calculations. We see that GGA leads to pla- 
nar structures, but LDA favors 3D structures for n > 7 
clusters. Thus, although the planarity of AU7 compared 
to the 3D structures of Ag7 or CU7 was attributed to rela- 
tivistic effects^, the observed planarity of Au clustersSiS 
is accounted for only using the GGA theory. We have also 
optimized the Au n structures in Figure ^ within GGA 
and LDA non relativistic pseudopotentials, resulting that 
the 3D structures become more stable energetically than 
the planar ones, except for n = 6 within GGA, as shown 
in the last row of Table [I] 



TABLE I: The binding energy difference AE b = E b (Au n ) - 
£i(Au*,), in eV, between the 2D and 3D isomers of the 
gold clusters represented in Fig. Q optimized within LDA 
and GGA xc-functionals using relativistic and non-relativistic 
(NR) pseudopotentials. 





Au 6 


Au 7 


Au 8 


Aug 




LDA GGA 


LDA GGA 


LDA GGA 


LDA GGA 




0.14 0.15 


-0.02 0.03 


-0.06 0.01 


-0.04 0.01 




-0.04 0.02 


-1.42 -1.11 


-1.64 -1.37 


-1.29 -0.98 



Table [n] gives the various energy differences (total, ki- 
netic, Coulomb, and exchange-correlation) between the 
second and hrst energy isomers of scalar relativistic Au 
and Cu clusters with 6 < n < 9 atoms. We will denote 
compact 3D structures with an asterisk, while its absence 
indicates planar and cage-like geometries. The geometry 
of the Au isomers is given in Figure Q The geometries 
of Cu6 and Cug are similar to those of Au6 and Aug re- 
spectively. For Cu„ with n = 7, 8, 9, the geometries 
of the two lowest energy isomers are both 3D, and are 
taken from our previous worki We can see in Table [H] 
that planar structures have smaller kinetic energy than 
3D isomers, and larger for LDA than for GGA. Adding 
kinetic and Coulomb energies, the 2D structures became 
more stable energetically than the 3D ones. On the other 
hand, the xc-energy is more negative (it contributes more 
to the binding energy) for 3D than for planar structures. 
Both effects are stronger within the LDA but, in the bal- 
ance of total energy difference, the loss of kinetic energy 
in planar structures dominates over the increase of xc- 
energy when using GGA, but the opposite occurs with 
LDA. Thus, as a whole, GGA (LDA) favors 2D (3D) 
structures of gold clusters. On the other hand, for Cu6 
and Cug isomers, which have the same geometry as Aue 
and Aug, the change in kinetic and Coulomb energies is 
not so noticeable as in gold. 

For Cu„ with n — 7-9, whose first and second isomer 
geometries are all 3D, the change in exchange-correlation 
energy, A xc , can be positive or negative. For Cu§, whose 
electronic properties can be described approximately by 
the spherical jellium model^, we see that the sum of 
changes in kinetic and Coulomb energies roughly can- 
cel, resulting that the change in the total energy is ruled 
by the change in xc-energy, as in the jellium model^. 

The loss of kinetic energy in planar gold clusters with 
respect to their 3D isomers was attributed to electron 
derealization^, but is not easy to reconcile that der- 
ealization with the simultaneous confinement in two di- 
mensions. On the other hand, the xc energy becomes 
less negative for planar configurations, which is also not 
clearly related to derealization and confinement in 2D 
gold clusters. We observe that the calculated average 
bond length, d av , is larger for 3D than for 2D isomers. 
Specifically, the difference between d av of 3D and 2D Au„ 
isomers with n = 6, 7, 8, 9, is (in A) 0.04, 0.11, 0.13, 0.09 
(0.05, 0.11, 0.09, 0.09) for the LDA (GGA) calculation. 

The kinetic energy of the electron gas in two dimen- 
sions is about one half of the 3D case at the same den- 
sity parameter, r s , while the exchange energy in 2D is 
slightly larger than in 3DiS^i. However, the correlation 
energy, at least in the RPA, is much larger in 2D than in 
cj-q46 Thjg consideration points again to the importance 
of good exchange-correlation functionals in dealing with 
the DFT structural description of gold clusters. 
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TABLE II: Total, kinetic, Coulomb, and exchange- correla- 
tion energy differences (AEi = £7j(Au£) — Ei(Au n )), in eV, 
between the second isomer (3D, Au^) and first isomer (2D, 
Au n ) of gold clusters with the structures of FigQ For Cu„ 
the first and second isomers are 3D, except for Cu6, where the 
two isomers are similar to Au6 and Aug. The Cu n geometries 
are taken from ref. 0- 





AE 


tot 


AE 


k i n 


AE C 


O u I 


AExc 




LDA 


CCA 


LDA 


GGA 


LDA 


GGA 


LDA 


GGA 


Au 6 


0.86 


0.88 


11.01 


10.63 


-9.62 


-9.42 


-0.53 


-0.33 


Au 7 


-0.15 


0.24 


12.49 


10.34 


-11.12 


-9.52 


-1.53 


-0.57 


Au 8 


-0.46 


0.04 


7.37 


7.03 


-6.53 


-6.80 


-1.30 


-0.19 


Aug 


-0.34 


0.07 


5.07 


4.63 


-4.50 


-4.52 


-0.91 


-0.04 


Cu 6 




0.03 




1.03 




-078 




-0.22 


Cu 7 




0.34 




0.20 




-0.98 




1.13 


Cu 8 




0.89 




2.10 




-1.98 




0.77 


Cu 9 




0.25 




-1.60 




1.87 




-0.02 



B. Magic cage-like structures of gold clusters 

Stable cage-like structures of gold clusters have been 
predicted recently for Au ^ 14 ^ , Au 2 g^, Au42^, and 
others 16 ' 17 . On the other hand, although all atoms of 
AU20 are at the surface, this cluster can be considered as 
a small piece of bulk fee gold^i2i±£. 

In this work we construct cage-like atomic structures 
starting with the Platonic solids with triangular faces 
— tetrahedron (Td), octahedron (Oh), and icosahedron 
(Ih) — which are those allowing compact planar pack- 
ing. By adding atoms at the intersections of fee planes 
in the triangles, we obtain the following sequences for 
the number of atoms: n(Td)=4 + 2m(m+2), n(0/ l )=6 
+ 4ra(m+2), and n(I^)=12 + 10m(m+2), where m = 
0, 1, 2,. . . is the number of atoms inserted in each edge. 
When we add a central atom to each new triangle, we 
obtain cage-like structures with n(T^) = 2+6(m+l) 2 , 
n(O h ) = 2+12(m+l) 2 , and 71(1/0 = 2+30(n+l) 2 . On 
the other hand, from the electronic point of view, a clus- 
ter with nearly free valence electrons is magic when it 
has filled electronic shells, having well defined angular 
momentum, as in the jellium model. As we look for 
an empty-cage with an approximately spherical surface, 
orbitals with radial nodes have to be excluded and the 
only allowed electronic shells are Is, lp, Id, If, . . . This 
leads to a magic number of electrons n e = 2(1 + l) 2 , 
where I is an integer numbei^. Assuming that each no- 
ble metal atom contributes with one valence electron, 
the equality n — 2(1 + l) 2 leads to the following magic 
neutral cage-like Au„ clusters (containing less than 1000 
atoms): n = 8 (T d ), 18 (O h ), 32 (I h ), 50 (O h )> 98(T d ), 
162 (Ih), and 578 (O^). Double anionic clusters should 
obey n + 2 = 2(1 + l) 2 and they appear at n = 6 and 
198 (Oh) atoms. Double cationic magic clusters are the 
solutions oi n — 2 — 2(1 + l) 2 . However, we exclude the 
Td clusters because they are far from spherical. Also 
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FIG. 2: (Color online) Cage-like (Au„) and space filling (Au^) 
equilibrium isomeric structures of neutral gold clusters with n 
= 18, 20, 32, 50, and 162, except Auig2 which is not optimized 
in this paper. 

AU20 can not be properly considered cage-like because 
it contains many internal bonds, as similarly occurs to 
smaller clusters, like Auis- Nevertheless, in the following 
we will test the estructural and electronic properties of 
Au„ clusters with 18, 20, 32, 50, and 162 atoms. 



TABLE III: Binding energy difference between the cage-like 
and compact structures of Fig|5| AEi = E B (A\l rl ) — Et(Au^), 
optimized with the LDA and GGA approximations. 



Auis 


AU20 


AU32 


Au 5 o 


LDA GGA 


LDA GGA 


LDA GGA 


LDA GGA 


AE b -0.07 0.01 


0.01 0.03 


-0.06 0.01 


-0.14 0.01 



We performed full relaxations of the initial cage-like 
magic structures and several compact geometries ob- 
tained by forcing initially some surface atoms inside 
these clusters. Fig [21 shows the equilibrium geometries 
with cage-like and compact structures, obtained after a 
non exhaustive search, and optimized with forces < 0.01 
eV/A) at the GGA and LDA levels. The cage-like equi- 
librium structures were proven to be meta-stable after 
performing an ab-initio molecular dynamics run at tem- 
perature of 600 K during 1000 steps, each of 2 fs. The 
binding energy difference between cage-like and compact 
structures is tabulated in Table ITTll We see that cage-like 
GGA structures are slightly more bound than the com- 
pact ones. Instead, LDA leads to compact structures 
without symmetry, except for Au2o- The true ground 
state is not known, however, and improved functionals 
could lead to compact, ordered or disordered, structures. 

As a further test, we calculate the relative stability of 
cage-like and compact isomers after loss or gain of one 
electron. Fig[3]shows the total energy difference per atom 
between cage-like and compact equilibrium structures of 
cationic, neutral, and anionic clusters with n=18, 20, 32 
atoms. We see that only the cation AujJ~ 2 and the anion 
Augo are still cage-like. Interestingly, the lowest energy 
isomer of both ionic clusters, Au^" 8 and Au^g, is not cage- 
like, contrary to the neutral Auis. 

In Fig 0] we compare the density of states (DOS) of 
the cage-like and amorphous structures of neutral Au n 



(n=18, 20). We see that the cage-like geometries show a 
more structured DOS, with well defined peaks. This is 
probably due to a higher geometrical symmetry (O^ for 
Auis and Td for AU20), as shown in a previous work for 
AusgiS. This fact is independent of the GGA or LDA 
level of theory. The LDA DOS profile tends to be shifted 
(to lower energies for AU20) with respect to the GGA 
one, and the HOMO-LUMO gap is smaller for LDA than 
for GGA. On the other hand the HOMO-LUMO gap is 
considerably larger for cage-like AU20 than for cage-like 
Auis, indicating that the former is much more stable 
than the later. 
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FIG. 3: (Color online) Total GGA energy difference per atom 
(in eV) between the lower energy isomer of compact and cage- 
like structures for cationic, neutral, and anionic clusters of 
gold with n=18 (stars), 20 (crosses), and 32 (circles) atoms. 



FIG. 4: (Color online) Density of states (DOS) of the lower 
energy cage-like and amorphous isomers of Au clusters with 
18 and 20 atoms, calculated at the LDA and GGA levels of 
theory. 



C. Static dipole polarizability 

The minimum polarizability principle states that any 
system evolves naturally towards a state of minimum po- 
larizability, but exceptions have been reported^. In Ta- 
ble llVl are given the results of our scalar-relativistic GGA 
calculations for the mean polarizability per atom of the 
two lowest energy states of the Au clusters reported in 
the subsections above. 

The static dipole polarizability of atomic Au comes 
to 20.53 a.u. (Bohr 3 ) for the scalar-relativistic calcula- 
tion, and 33.06 a.u. for the non-relativistic calculation, 
which is a clear manifestation of the relativistic size con- 
traction of gold 2 *. Our Au polarizability is smaller than 
other calculated and experimental values quoted in the 
literature2i£2iSi£i. Different estimates of the experimen- 
tal value are 30±4 a.ui£ and 39.1±9.8 a.u. 52 . A quan- 
tum chemical CCSD(T) calculation gives 36.06 a.ui^, 
and a recent CASSCF/CASPT2 relativistic calculation 
gives 27.9 a.uS. 

For noble metal dimers, recent relativistic and non- 
relativistic calculations^ 5 - found that the reduction of the 
polarizability, due to the relativistic contraction effect, 
amounts to 39.8 % for Au 2 , 15.8 % for Ag 2 , and 6% for 
Cu2 . From the TDDFT study of Castro and coworkers^, 
the static polarizability of Au„ clusters up to n = 4 are 



affected up to ~ 2% by the inclusion of the spin-orbit 
term. 

From the TDDFT calculations of Castro 5 ^, the static 
polarizability of the 3D tetrahedral isomer of Au4 is 2.5% 
higher than the one for the planar T>2h ground state. We 
obtain a larger difference (15.5 %) between our two pla- 
nar Au4 isomers (a rhombus and a triangle with another 
Au on top, see refi. Our result for the Au„ polarizability 
with n = 2-9, is similar to the one calculated by Zhao and 
coworkers 5 ^ using a finite field perturbation method like 
ours. However, the polarizability per atom obtained by 
these authors for planar Auis and Au 2 o is much higher 
than ours for cage-like and compact isomers. On the 
other hand, the difference between GGA and LDA polar- 
izabilitics of planar Au„ clusters (n = 2-20) calculated in 
ref. 03 is less than 2%. 

We see in Table HVl that the average polarizability per 
atom of Au n clusters is remarkably constant, with small 
odd-even effects for n < 7. Up to n — 6 the first and 
second isomers are planar, and the larger polarizability 
of the second isomer is due to the larger average Au- 
Au bond length. The polarizability of the Au£ 3D iso- 
mer is smaller than that of the planar isomer, due to the 
compact pentagonal bipyramid geometry, and it consti- 
tutes one of the exceptions to the minimum polarizability 
principle 51 . The cage-like structures seem to be another 
exception to that rule. 

The average mean polarizability per atom for the first 
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TABLE IV: The third and fourth columns give, respectively, 
the mean polarizability per atom a/n of the first and second 
isomers of Au n clusters. The last column gives their differ- 
ence. For comparison, the first column gives the a/n value 
calculated by Zhao and coworkers for planar estructures— . 
All polarizabilities are in atomic units (bohr 3 ). 



valence, n = number of atoms) as 







a/n 


a/n 




n 


Ref. 57 


1st isomer 


2nd isomer 


Aa (%) 


2 


29.4 


28.3 


29.0 


2.5 


•i 
o 


or . i 






18 9 


4 


32.7 


31.5 


37.2 


17.9 


5 


34.2 


32.9 


46.6 


41.6 


6 


34.2 


32.7 


40.0 


22.5 


7 


34.3 


33.9 


30.9 


-9.0 


8 


36.2 


34.7 


36.1 


4.2 


9 


38.5 


37.8 


48.1 


27.3 


18 


42.7 


32.1 


31.4 


-2.4 


20 


43.7 


35.2 


32.9 


-6.4 


32 




33.9 


33.5 


-1.3 



isomer of Au„ clusters in Table lIVl is a av /n — 33.57 a.u.. 
The jellium model for atomic clusters of monovalent met- 
als fulfills the relation ctj e i — r 3 n, where r s is the radius 
per electron (in a.u.) of the bulk metal. For Au, with r s 
= 3.01 and the GGA, we obtain the ratio a = a av /aj e i 
— 1.23. We can calculate a similar ratio for Li, Na, Cu, 
and Ag clusters using results from experiments and other 
calculations. Using the experimental valuesSiS for the 
dipole polarizability of Li (r s = 3.25), Na (r s = 3.93), and 
Cu (r s =2.67) clusters with n < 10, we obtain the ratios 
a exp (Li) = 2.26, cv e:E p(Na) = 2.11, and a exp (Li) = 2.0 
for Li, Na, and Cu, respectively. From the GGA polariz- 
ability values calculated in ref. |H3 for Li and Na clusters 
with n < 10, we find acGAiX^) = 2.47 and o.gga (Na) = 
1.81. For Cu clusters in the same range of sizes, using the 
different GGA and LDA polarizability values reported 
by Yang and JacksonSi, one has Sgg^Cu) = 1.93-1.98, 
and 0(lda(Cu) — 1.79. Using recent calculations-^ for 
the polarizability of Ag clusters with 2 < n < 8, results 
in acG^Ag) = 1.69. As a whole, we see that the ratio 
aav/ajei is considerably smaller for Au than for Li, Na, 
Cu, and Ag clusters. 

The enhancement of dipole polarizabilty over the clas- 
sical jellium model (Mie value <x, e z) is directly propor- 
tional to the fraction of electronic charge that extends 
beyond the positive background in the field free system 
(spill-out) 1 . The smaller spill-out for gold with respect 
to silver and copper can be attributed tentatively to the 
relativistic contraction of the electronic cloudSiSi. 

We can explore a little more our calculated polarizabil- 
ity values using an extended Thomas-Fermi- Weizsacker 
(TFAW) jellium model 6 ^, which predicts the mean po- 
larizability of a cluster with n e = vn valence electrons {v 



a jel[ 



1/3 



(1 



^%V3 }] (2) 



In this expression, d(r s ) is the image plane position (the 
centroid of the induced electron density for the flat metal 
surface), and ot\ and a are two coefficients dependent of 
the parameter A, which take into account the weight of 
inhomogeneity-density correction (Weizsacker term) to 
the Thomas-Fermi kinetic energy. In the interval 3 < r s 
< 4, it results cvi(r s ) ~ -0.1 a.u. independently of the 
A value (see Fig. 1 of ref. :62), but a (r s ) is strongly 
dependent of A. In terms of a reduced polarizability, 
defined by a re d = (a/a je ; — l)n. 2 / 3 , we can write eq. 
as 



a red = An 1/3 + AB 



where A 



3 d(r s ) 
7Aj3 r , i 



and B 



(3) 

^T^TTTJ- % fittin s 

our data in Table IIVI to eqn. [21 we obtain a re d = 
1.1571 1 / 3 — 1.26, and using d(Au)=1.34 a.u. for the image 
plane distance of gold 6 —, there results an effective va- 
lence v = 1.56, and a ratio ai ^ s \ — -1.46. For ai(Au) ~ 
-0.1 we obtain an( Au) = 0.07. Extrapolating the results 
in Fig. 1 of ref. such a small ao(Au) value corre- 
sponds to an extremely small value of the parameter A, 
which means that the contribution of inhomogeneity cor- 
rections to the kinetic energy is very small. This result 
agrees with our conclusion in subsection IIII Al about the 
delocalization of valence electrons in planar gold clusters. 
On the other hand, the effective valence v = 1.56 reflects 
in some way the screening of the dipole response due to 
the <i-electrons£&. 

We test also an empirical linear relation between the 
cubic root of the mean polarizability per atom, a 1 / 3 /n, 
and the inverse of the ionization potential per atom, 
lp 1 /n. That relation was tested in ref. |6(3 for Li and 
Na clusters with 2-10 atoms, resulting for both cases (us- 
ing calculated GGA values, in a.u.) in a proportionality 
constant close to unity, with a linear correlation coeffi- 
cient 0.995. For the lower energy isomers of Au clusters 
reported in Table ITVI there results a proportionality con- 
stant 0.925 and a correlation coefficient of 0.984. 

As probed by Yang and Jackson 61 , temperature effects 
are a possible source of discrepancy between calculated 
and measured polarizabilities, because calculations are 
carried out for K while experiments are conducted at 
finite temperatures. The existence of a permanent elec- 
tric dipole in a cluster adds the following temperature- 
dependent term to the effective polarizability 64 : 



a eff 



J£_ 

3fcT' 



(4) 



where fi is the dipole moment and k is the Boltzman con- 
stant. The dipole contribution is important at low tem- 
peratures for clusters with permanent dipole moments. 
Using the dipole moments calculated for our GGA lower 
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energy isomers, the correction of eqn. 0] to the mean 
polarizability per atom is (in a.u) 1.39, 0.04, 0.74, 2.33, 
and 0.09, for clusters with 3, 4, 5, 6, and 20 atoms, re- 
spectively. This small correction still allows one to dis- 
criminate the polarizability of planar and cage-like gold 
clusters from their isomers. 



IV. CONCLUSIONS 

We obtain that, using non-relativistic pseudopoten- 
tials, both GGA and LDA predict the onset of three di- 
mensional cluster structures already at n — 6 for Cu„ 
and Ag„, and at n — 7 (6) for Au„. This result changes 
by considering scalar-relativistic pseudopotentials within 
GGA, resulting in planar Au„ structures up to n = 11. 

From our scalar-relativistic results for the two lowest 
energy isomers of Au„ and Cu„ with n = 6-9 atoms, we 
find that planar structures have smaller kinetic energy 
than 3D isomers, and this effect is much larger for gold 
than for copper clusters. Adding kinetic and Coulomb 
energies, the 2D structures became more stable than the 
3D ones, and this effect is more noticeable for LDA than 
for GGA. On the other hand, the xc-energy is more nega- 
tive (contribute more to the binding energy) for 3D than 
for planar structures, and this effect is notably enhanced 
within LDA. Thus, in the total energy balance, kinetic 
energy loss favors planar GGA structures, but xc-energy 
favors LDA 3D structures. 

As a second step, we constructed clusters having only 
surface atoms, and with O^, T^, and 1^ symmetry. From 
those, assuming one valence electron per atom, we select 
the ones having 2(l+l) 2 electrons, which correspond to 
the filling of a spherical electronic shell formed by node- 
less one electron wave functions. We obtain, by means 
of scalar-relativistic GGA calculations, that these cage- 



like structures for neutral Auis, AU20, AU32, AU50, and 
Aui62 are meta-stable after moderate (600 K) constant 
temperature molecular dynamics. However, after addi- 
tion or substraction of an electron, only the anion Au^q 
and the cation Au^~ 2 remain cage-like. 

Finally, we calculate the static polarizability of the two 
lowest energy isomers of Au„ clusters, which are planar 
(n = 2-9) and cage-like (n = 18, 20, 32) for the first 
isomer and planar (n = 2-6) or space filling 3D (n = 7- 
9, 18, 20, 32) for the second isomer. In the range n = 
2-9, the polarizability per atom is smaller for the first 
isomer than for the second, with the exception of n = 7, 
confirming the empirical rule of minimum polarizability. 
The contrary occurs for cage-like structures, with larger 
a than their space-filing isomers. 

We fitted the polarizability of the first isomer of these 
gold clusters to a semi-empirical relation between the 
cluster dipole polarizability and its size, which involves 
the effective atomic valence and the kinetic energy due 
to homogeneous and to inhomogcneous density compo- 
nents. From that fit we extract a very small value for the 
kinetic energy component due to inhomogeneous density, 
which suggest a delocalized character of the valence elec- 
trons involved in the dipole response. We also obtain 
an effective valence charge of 1.56 electrons, reflecting a 
dipole response with strong screening of the d-electrons, 
as already reported previously^. 
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